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We introduce a one-parameter deformation for one-dimensional quantum lattice models, 
the hyperbolic deformation, where the scale of the local energy is proportional to cosh Xj at 
the jih site. From the quantum-classical correspondence with a two-dimensional classical 
system, it is expected that the deformation does not modify the ground state conspicuously, 
if there is a finite excitation gap. Under this situation, the excited quasi particle is weakly 
confined around the center of the system. We derive scaling relations among the mean-square 
width {w 2 } of confinement, energy correction to the excitation gap A, and deformation 
parameter A. This scaling is useful for the extraction of the bulk excitation gap in the limit 
A = 0. 

§1. Introduction 

The ground state of an infinitely large quantum system that has a nonzero exci- 
tation gap is quite different from those of gapless systems in its correlation properties. 
Precise estimation of the excitation gap is therefore one of the theoretical interests 
in condensed matter physics. For one-dimensional (ID) quantum systems, the den- 
sity matrix renormalization group (DMRG) method^ 'I 3 -* is often used for the 
gap estimation, since it gives low-lying eigenvalues with high numerical precision. 
Since the target systems of DMRG are chiefly finite-size ones under open boundary 
condition, finite-size scaling (FSS^P 1 '® has been employed for the subtraction of the 
finite-size effect, which includes boundary corrections. 

Energies of the ground state and the first excited state are necessary for the de- 
termination of the excitation gap. The ground-state energy E ® for a finite system 
of size L — typically up to the order of L ~ 1000 — can be calculated numer- 
ically. Under the open boundary condition, an appropriate tuning of the system 
boundary is important for a rapid convergence of E^ /L with respect to L, where 

= ]hnjj_ ¥oo Ejp / L is the ground-state energy per site in the thermodynamic 
limit. To speak from a numerical viewpoint, there are several ways to obtain a faster 
convergence of with respect to L. Since the number of bonds of the L-site system 
is L — 1, E^ /(L — 1) could be better than E^ JL. To observe the bond energy at 
the center of the system or to calculate (E^) 2 — E^)/2 is more efficient to avoid 
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the boundary energy corrections. Further improvements can be achieved by a fine 
tuning of interactions near the system boundary. For example, the introduction of 
weak bonds near the boundary^® realizes a very fast decay of the boundary effect. 
To introduce sinusoidal modulation to the local energy scale is also useful for the 
suppression of the boundary effect P 

Precise estimation of the excitation gap A is more difficult than that of e^°\ 
even when the quasi-particle picture naturally holds for the elementary excitation. 
This is because the extrapolation process 

A= lim (E^-EfA (1-1) 

L— s-oo \ / 

requires energy of the first excited state E^ , which is easily affected by the reflec- 
tion of the excited quasi particle at the system boundary. Such an effect may cause 
nonnegligible L-dependence on Ejp — E^ . Thus, it is important to tune the inter- 
actions near the system boundary, so that the kinetic energy of the excited quasi 
particle rapidly (or regularly) converges to zero with respect to L.ED>in)ilI3 Gener- 
ally speaking, it is not easy to perform the boundary tuning, which works efficiently 
for both ground and excited states. Obviously, the introduction of very weak bonds 
near the system boundary creates a spurious localized excitation near the boundary. 

In this article, we consider a way of suppressing the boundary reflection effect, 
by weakly confining the excited quasi particle around the center of the system, pre- 
serving uniformity in the ground state as much as possible. To clarify the situation, 
let us consider the Hamiltonian 



on ID lattice labeled by the site index j, where /t - J+1 and g - denote neighboring in- 
teraction and on-site term, respectively. We introduce a one-parameter deformation, 
the hyperbolic deformation, to the above HamiltonianP^} The deformed Hamiltonian 
is written as 

H(X) = ^2 cosh A -? hjj+i + cosh A ( J ~ I) & ' ( 1-3 ) 

j j 

where A is the deformation parameter. The ground state of this deformed sys- 
tem is expected to be nearly uniform, although the interaction strength is position- 
dependent. There is a geometrical reason for this uniformity, which we explain in 
the next section. 

Since the energy scale increases with respect to \j\ in Eq. (1-3), the excited quasi 
particle cannot reach the system boundary when XL is sufficiently large. This con- 
finement gives a finite size correction to the energy of the first excited state E^ (A), 

where the correction E^ (A) — E$ (0) comes from the uncertainty relation and its 
analytic form can be estimated. In order to detect this correction quantitatively 
for various gapped systems, we propose an effective one-particle Hamiltonian for 
the excited quasi particle. As we show in §3, there is a scaling relation among 
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L, E£>(\) - E^'(X), and A under this quasi-particle picture. This scaling can be 
applicable for those systems that possess a nice quasi-particle picture, such as the 
Hubbard modeP^^ and inte ger-spin chains. Conclusions are summarized in the 
last section, and the applicability to the S = 1 antiferromagnetic Heisenberg spin 
chain is discussed. 

§2. Uniformity under the hyperbolic deformation 

The real- or imaginary-time evolution of a ID quantum system is related to a 
2D classical system through the so-called quantum-classical correspondence 113 .CHI 1 'El 1 
Let us consider the correspondence for a case where a classical system is on a curved 
2D space. We focus on the hyperbolic plane, which is a 2D space with a constant 
negative curvature. Suppose that there is a uniform action of the classical field on this 
curved surface. What does the corresponding ID quantum Hamiltonian look like? 
Let us consider this problem from the viewpoint of the imaginary-time evolution. 



/ s 
/ ^ 
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Fig. 1. Imaginary-time axis shown by the vertical line and horizontal equal-time curves on the 
hyperbolic plane drawn inside the Poincare disc. 

Figure 1 shows the 2D hyperbolic plane drawn inside the Poincare disc. All 
the geodesies are represented by arcs, which are perpendicular to the border circle, 
including straight lines that pass through the center of the disc. Let us choose 
the vertical line as the imaginary-time axis, where we denote the time by r. All 
the geodesies that are perpendicular to this imaginary-time axis can be regarded as 
equal-time curves. Consider a quantum state |^(t)) on an equal-time line. If the 
classical action in the lower half of the hyperbolic plane is uniform, and if there is 
no symmetry breaking such as dimerization, the state |^(r)) is also translationally 
invariant. 

Let us consider a tiny evolution from r to r + St 

\$(t + St)) = U[6t] |#(t)) , (2-1) 

where W[(5r] represents the imaginary-time evolution. Although we have assumed 
that both ^(t)) and \&(r + St)) are translationally invariant, U[5t] is not. This 
fact can be understood geometrically by considering the distance between two points 
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(x, t) and (x, r + 8t) on the hyperbolic plane, where x represents the spatial coor- 
dinate. The distance is an increasing function of \x\ as shown in Fig. 1, and can 
be written as (coshvx) 5t, where v is a function of the scalar curvature of the 
hyperbolic plane. If it is possible to represent lA\br\ in the form exp(— brH) , the 
corresponding Hamiltonian H is also position-dependent, and might be written with 
an integral 

H = j (cosh.vx)h(x)dx (2-2) 

of a local operator h(x). This is a rough sketch of the hyperbolic deformation of a 
quantum Hamiltonian in the continuous ID space. 

Let us consider a construction of the hyperbolic deformation in the discrete 
space. We introduce lattice points at x = aj, where a is the lattice constant, and j 
the site index. (See Appendix A for the correspondence for the one-particle case.) 
We also define a new dimensionless parameter A = va\ we have 

cosh vx = cosh vaj = cosh Aj . (2 - 3) 

A discrete analogue of H in Eq. (2-2) is then given by Eq. (1-3), which is H(X) = 
Y^j cosh Aj hj • +1 + Ylj cosh A (j — |) <7 - . When A = 0, the Hamiltonian H(X = 0) 
coincides with the translationally invariant one in Eq. (IT). From the construction of 
.ff (A), it is expected that the ground state of this deformed system is nearly uniform, 
although the interaction strength in H (A) is position-dependent. This assumption of 
the uniformity can be explicitly verified for classical statistical systems on hyperbolic 
lattices!^* Moreover, it was observed that the bond energy of the ground state of the 
5 = 1/2 Heisenberg chain has very small position-dependence when the hyperbolic 
deformation is imposed.^ 1 

§3. Finite-size scaling for the excitation energy 

Compared with the ground state, the effect of the hyperbolic deformation on the 
excited states is nontrivial. This is because quasi particles are subjected to the local 
interaction with the energy scale proportional to cosh Aj. We analyze the behavior of 
the excited quasi particle using an effective one-particle model. Under the hyperbolic 
deformation, we can naturally expect that the effective Hamiltonian for the quasi 
particle also takes the form of the hyperbolic deformation 

H F (X) = ~t^2 cosh A -? ( 4s+! + 4+iS') ~ coshX ~ \) > I 3 " 1 ) 

j j 

where and Cj represent, respectively, creation and annihilation of the quasi particle. 
They are either Fermionic or Bosonic, since we consider the zero-particle vacuum |0) 
and one-particle states !<?) only, and therefore, the statistics of these operators does 
not matter in the following analysis. The letter 'F' in H F (X) denotes that H F (X) is 

*' The distance (cosh vx) St diverges in the limit \x\ — >• oo. This divergence can be controlled 
by considering finite width R and keeping \x\ < R. 
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quadratic, and is free from interaction terms. For the moment, we assume that the 
system size is sufficiently large. 

We have denoted the hopping amplitude of the quasi particle by t. When the 
Hamiltonian is undeformed (A = 0), the parameter fi can be interpreted as the 
chemical potential. Since we have introduced H F (X) as an effective Hamiltonian 
for the gapped system, we consider the case when \x < — 2t cosh(A/2). Under this 
condition, the ground state is trivially the zero-particle vacuum |0) with zero energy 
E^\\) = 0, where L is the system size as was introduced in §1, and the lowest 
excited state is the one-particle state \&) with positive energy. It should be noted 
that we have assumed the translational invariance of the ground state of H(X) in 
Eq. (1-3), when we extract the effective Hamiltonian H F (X). 

Even for the toy model in Eq. (3-1), the distribution of the quasi particle for 
the lowest-energy one-particle state is nottrivial. The second term in Eq. (3-1) 
is relevant to the confinement effect. Introducing the wave function \Pj = = 
(0|c J -|$ r ), we have 

4 1 ' (A) &j = -t cosh Xj - t cosh A (j - l) - n cosh A (j - |) ^ , (3-2) 

which is nothing but the one-particle Schrodinger equation. (See Appendix A con- 
cerning the continuous limit.) Note that the single-particle excitation energy when 
L = oo can be written as 

ZUA) = e£\X) - E®{\) = £«(A) , (3-3) 

which has the small-A limit 

A=limA oc (X) = -2t- f i. (3-4) 

A->0 

Our aim is to find an efficient finite-size scaling, which captures this limit correctly. 
It is easy to solve Eq. (|3-2p numerically when the system size L is finite, up to L of 
the order of 10 2 ~ 10 4 . Thus, what we have to do is to take the double limit 



A = lim 



4 1) (A)-4 0) (A) = lim A L (X) (3-5) 

A — Hj 



somehow by means of a scaling plot. If there is an efficient finite-size scaling for 
Al(\), it could be useful for gap estimation of many-body models by means of the 
hyperbolic deformation in Eq. (1-3). As shown in the following, the double limit 
A — > and L — > oo can be stably taken by keeping (L + 1) 2 A constant. 

Hereafter, we choose t as the unit of energy, i.e., t = 1. From the relation in 
Eq. (|3-4p . the parameter \i can be represented as/i = — 2t — A = — 2 — A. In the 
following, we specify A instead of /x to simplify the formulations in the finite-size 
scaling. Figure 2 shows the occupancy pj = ^j^j of the lowest-energy one-particle 
state when A = 1 and L = 100. We have labeled the sites with j ranging from 
j = -L/2 + 1 to j = L/2. When A = 0.01 or 0.02, the Gaussian-like distribution 
is observed around the center of the system j = 1/2. The distribution can be 
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Fig. 2. Site occupancy pj = tyj&j under the conditions L = 100, t = 1, and /i 
in Eq. (3-4) is equal to unity. 



-3. The gap A 



characterized using the mean-square width 



L/2 , , 2 

<- 2 >= £ (J-I) Pr ( 3 - 6 ) 



j=-L/2+l 



Consider the case where the parameter A is relatively small. As we have exam- 
ined in Fig. 2, confinement of the quasi particle appears when L is sufficiently large 
and XL S> 1 is satisfied. Under this condition, we can estimate the A dependence 
of (w 2 ) from the uncertainty relation. The typical energy to put a localized wave 
packet around the jth site is of the order of 



A cosh A (j - \) ~ A 



(3-7) 



where we have truncated the Taylor series by the second order in j — ^. Replacing 
(j — i) 2 by (w 2 ), we can roughly estimate the potential-like energy 

U~A + A — (w 2 ) (3-8) 

for the weakly confined state. From the uncertainty relation between position and 
momentum, the kinetic- like energy T = p 2 /2m can be estimated as 

(3-9) 



2 (w 2 ) ' 



where b is a constant of the order of t. Minimizing the sum E = U + T, we obtain 
the relation 

( w 2\ „ (3 . 10) 

X 1 VAX 



after a short calculation. Similarly, we can express the energy correction as E 
VbAX, which is caused by the confinement. 
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On the other hand, when A = 0, the confinement effect disappears. In this case, 
the quasi-particle wave function is exactly given by 



3 V L + 1 



2 j 
cos ■ 



L + l 



7T . 



(3-H) 



which represents a free particle inside the area of width L + l. Therefore, (w 2 ) is of 
the order of (L + l) 2 , where the proportional constant c is calculated as 



lim 



L->oc (L + l) 2 

The energy ^^(0) = E^\o) is obtained as 



f 1 x 2 cos 2 7rx dx = y ^ . (3-12) 

7-1/2 12 



-2t cos 



7T 



L+l 



-2 cos 



7T 



L + l 



+ 2 + Z\, 



(3-13) 



and therefore, the finite size correction to the true gap A = lim L ^ 00 A L (Q) is derived 
as 



A L (0) - A 



IT 2 „ / 1 \ 7T 2 C 



(L+l) 2 







(3-14) 
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Fig. 3. Scaling plot for (w 2 ). The solid line and 



UL+l) 

dotted line, 



respectively, correspond to (w 2 )/(L + 



l) 2 = c and {w 2 )/(L + l) 2 = Vb/y, where c = (1 - 6/tt 2 )/12, & = 1/2, and y = VZ\(L + l) 2 . 

These observations suggest that (A L (X) — A) (w 2 ) is a slowly varying function 
of A, and that the scaling relation 



(L + l) 2 



/ \VA\{L + iy 



(3-15) 



is satisfied. The function f[y] satisfies /[0] = c and f[y] ~ Vb/y when y is sufficiently 
large. Figure 3 shows the relation between (u> 2 ) and VAX(L + l) 2 for A = 0.5 and 
A = 1. The plot supports the presence of the scaling function /[y]. In the same 
manner, it is expected that the energy correction A L (X) — A satisfies the scaling 
relation 

(L + l) 2 (A L (X) -A) = g [VA\{L + l) 2 ] , (3-16) 
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where g[0] = ir 2 and g[y] ~ Vby when y is sufficiently large. Figure 4 shows the 
relation between (L+ l) 2 (A L (X) — A) and y/~AA(L + l) 2 . Clearly, there is a scaling 
function g[y] also for A L (\) — A. In summary, an advantage of hyperbolic defor- 
mation is that the effect of the system boundary disappears when \/AA(L + l) 2 
is sufficiently large. This property would make the numerical determination of A 
easier. 




10 u 



10 
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Fig. 4. Scaling plot for A L {\) — A. The solid Tine and dotted line, respectively, mean (L + 
1) 2 (A L (\) -A) = n 2 and (L + 1) 2 (A L (\) - A) = Vby. 



§4. Conclusion and discussion 



We have introduced the hyperbolic deformation to ID quantum lattice models. 
The quantum-classical correspondence suggests that the ground state is not strongly 
modified by the effect of deformation, while the excited quasi particle is weakly 
confined around the center of the system. The mean-square width of the quasi- 
particle distribution is calculated for a particle in the ID lattice. As a result, we 
have obtained the scaling relation Eq. (3-16) for the finite-size correction to the 
excitation gap. 




10' 10' 10 J 

Fig. 5. Scaling plot for A L (\) - A S = 1 Heisenberg spin chain. 



As a confirmation of the scaling relation in realistic systems, we calculate the 
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excitation energy of the hyperbolically deformed 5 = 1 antiferromagnetic Heisenberg 
spin chain, where h^ +1 in Eq. (1-3) is given by Sj ■ S - +1 , using the DMRG method. 
Figure 5 shows the calculated result plotted in the same manner as in Fig. 4, where 
the best fit is obtained when A is 0.410485. Numerical details and further investiga- 
tion would be presented elsewhere. It is thus expected that the scaling plot proposed 
in this article is useful for the numerical estimation of the excitation gap of various 
ID correlated systems, such as the higher S spin chains. In the case that the excited 
quasi particle occupies several sites on the lattice, one has to introduce an additional 
parameter d and replace (L + l) 2 by (L + d) 2 to reduce higher-order corrections. 
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Appendix A 

Continuous limit of the one-particle Hamiltonian 



Let us check the continuum limit of Eq. (j3-2|) . Substituting the relations x = aj, 
A = ai>, and the correspondence 9^ = &(aj) = &(x), we obtain 

E \P(x) = — t cosh v(x — — ) cosh v— \&(x + a) + ^(x — a)] 

—t sinh u(x — — ) sinh u— \W{x + a) — ^(x — a)] — /x cosh v (x — — ) \P(x) 

(A-1) 

h 2 

after some algebra. Expressing the hopping amplitude as t = 2ma $ , chemical poten- 

h 2 v 2 

tial as \x = — U — -^-g — 2t, and taking the limit a — > 0, we obtain a differential 
equation 

EV(x) = 



h2 8 , 9 rr , ' 

— cosh v x — — V U cosh vx 

2m ox ox 



9(x) . (A-2) 



The first term in the r.h.s. is the deformed kinetic energy, and the second term is 
a kind of trapping potential when U > 0. The Lagrangian that draws the above 
equation from the stationary condition is 



d 



h 2 d&*d& rrr ,' 



2m dx dx 



(A-3) 



where we have introduced the letter t for the imaginary-time variable, and where 
we have used the unit that satisfies h = 1. Note that the timelike variable r in 
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Eq. (2.1) is related to t by the relation (cosh 1/2;) dt = dr, and in the x-t plane, the 
Lagrangian can be represented as 



C'(9*,d T 9*, d x ¥*,9,d T ¥,d x V) = coshra 



H — — — — + & 



dr 2m dx dx 



The action of the system is therefore written as 



(A-4) 



d 

<fr*—y]/ + h{x) 

OT 



where the local Hamiltonian h(x) can be simply represented as 

~ ti 2 d&*d& TTTifT 
h(x) = - — — + U&*& 



2m dx dx 



(cosh vx) drdx , 

(A-5) 

(A-6) 



under this coordinate. 
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